(N 



o 

CD 



Isotropic turbulence in the dark fluid universe with 
inhomogeneous equation of state 

M. J. S. Houndjo^"'^!], M. E. Rodrigues^^j], C. Ainamon ^'^^ J. B. Chabi Orou^''^'^ and R. 

Myrzakulovt^j] 

" Departamento de Engenharia e Ciencias Naturais - CEUNES - Universidade Federal do Espirito 

Santo - CEP 29933-415 - Sdo Mateus/ ES, Brazil 
^ Institut de Mathematiques et de Sciences Physiques (IMSP) 01 BP 613 Porto-Novo, Benin 



o 

04 ■ Universidade Federal do Espirito Santo - Centro de Ciencias Exatas - Departamento de Fisica 



Av. Fernando Ferrari s/n - Campus de Coiabeiras CEP29075-910 - Vitoria/ES, Brazil 
\ Universite d'Abomey Calavi - Faculte des Sciences et Techniques (FAST) 

, Departement de physique - 01 BP : 4521 Cotonou 01 

^ Eurasian International Center for Theoretical Physics 

L.N. Cumilyov Eurasian National University, Astana 010008, Kazakhstan 

Abstract 

We investigate the turbulence effect in dark ffuid universe with linear inhomogeneous equation of 
«j ■ state. Attention is attached to two physical situations. First, we perform the perturbative analysis of 

■ turbulence and check its effects around the Big Rip. Later, treating the turbulence energy density as 

a part of total dark fluid, we study the stability of the system. The result shows that the stability is 
04 ' achieving as the energy density of turbulence decreases, changing into heat (the radiation), in perfect 

_ agreement with the avoidance of the Big Rip. 
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a : 1 Introduction 

There is a lot of interest in the study of the nature of dark energy, responsible for the acceleration of the 
cosmic expansion IHIl], and which is expected to have strange properties, like negative pressure and/or 
negative entropy. The current acceleration of the universe is realized when the parameter of equation of 
state (EoS) is such that uj < —1/3. More precisely, —1 < u < —1/3 corresponds to quintessence region 
where the expansion of the universe is accelerated, a, but not super-accelerated H > 0. When uj < —I, 
the universe is in a phantom region, which usually ends by finite-time future singularities [5]-|23j. 
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Recently Brevik et al analyse the possibility of whether there can be a turbulence micro-structure 
superimposed on the dark fluid [24 . In such situation, in order to maintain the macroscopic isotropic, 
attention is attached to the isotropic version of turbulence. An important aspect of this issue is that 
turbulence implies a loss of kinetic energy, changing it into heat. At macroscopic level, this turbulence 
is often treated as bulk viscosity. So, at microscopic level, within turbulence consideration, the length 
scale associated with the macroscopic viscosity may be replaced by the Kolmogorov microlength scale, 
conventionally denoted by ij [24| . In that paper, authors assumed homogeneous equation of state and 
studied the influence of the turbulence on the evolution of the universe, mainly, in it dark sector. 

In this paper, the goal remains the same, of studying the effects of the turbulence, but assuming 
inhomogeneous equation of state. We assume the EoS p — Lo{t)p + A(t) in which a;(t) and A(i) depend 
linearly on time. We also consider in this paper an oscillating w(t). This particular king of EoS is one 
alternative amongst a variety of possibilities, proposed to cope with the general dark energy problem 
[2S1[55], and developed for various purposes in p7]-[57]. 

Following the same steps of [24] we include the turbulence effect by adding a constant fraction, denoted 
/, to the laminar ordinary energy density in the first equation of FRW. Then, we perform a perturbative 
analysis to the first order, assuming / to be a small quantity and the correction term to uj{t) and K{t) 
have been found. Other interesting feature of considering this inhomogeneous EoS is that it generalizes 
the Brevik's results [24], at least, when FRW equations with the fraction / of the energy as turbulence 
energy are assumed. 

Moreover, instead of only considering the above perturbative treatment, we also consider the total 
dark energy contribution as a laminar part plus a turbulence energy part. Then, we study the stability 
of the system. The existence of an attractor critical point is investigated which should correspond to 
the stability of the system. Within the assumption that the dynamics of the system essentially depends 
on the evolution of the turbulence and the radiation, we find that stability corresponds to a decreasing 
turbulence energy density and an increasing radiation energy density. 

The paper is organized as follows: in section 2, we present the Kolmogorov's isotropic turbulence 
within inhomogeneous EoS. Section 3 is devoted to the study of turbulence effects as a fraction of the 
energy in the FRW equations, and this, for various assumptions of parameters Lo{t) and A(t). We inves- 
tigate the stability of the system in the section 4 where turbulence is assumed as an additive component 
to the total dark fiuid. We present the conclusion in section 5. 
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2 Kolmogorov's isotropic turbulence within inhomogeneous equa- 
tion of state 



In general, turbulence leads to a loss of kinetic energy into heat. As put out in [24l [38l |40], one 
can replace the macroscopic bulk viscosity with a microscopic shear viscosity, yielding the Kolmogorov 
length, by 



where v is the kinematic microscopic shear viscosity and e is the dissipation per unit time and unit mass. 

Let us consider I the external scale of the turbulence, and therefore, l/l the corresponding wave 
number. We also assume that the large eddies move around with only a little dissipation of energy. It 
is known from the second hypothesis of Kolmogorov [241 138) that in an isotropic region the motion is 
entirely determined by friction and inertia. There occurs a continuous flux of energy transferred by means 
of a hierarchy of eddies corresponding to the dissipation e. Let A denote the size of an eddy, with k — l/X 
the corresponding wave number. The equilibrium range is that for which all memory of the flow is lost, 

fc>l//. (2) 

If ux characterizes the typical velocity of an eddy of size A, the internal Reynolds number is written as 
ReA ^ \u\/v. For increasing values of fc, ReA decreases. When ReA ^ 1, dissipation becomes important, 
corresponding to the condition which implies the Kolmogorov length ([1]). For high value of Reynolds 
number, the wave numbers \/l and \/rj are widely separated, and there exists an inertial subrange 
characterized by [21] 

7«A:«-, (3) 

where the fluid behaves like a non-viscous fluid. In the inertial subrange the formula of the spectral 
energy density, conventionally called E(k), is 

with a w 1.5 being the Kolmogorov constant. 

For practical purposes the von Karman interpolation formula for Eik), linking the region of small k to 
the region of high k, is useful in order to calculate the total energy density E by integrating over all wave 
numbers (see Ref. [39j ) . In this paper, we abstain from going into further detail about this issue, and just 
focus our attention to the empirical decay law of the isotropic turbulence. Based on grid experiments in 



wind and water tunnels, it turns out that the mean kinetic energy iu^(t) decays as 

^[2(i)cxr6/^ (5) 
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see Refs. [HI |42], as well as the theoretical treatment in ^39j. This relationship shall be used in what 
follows. 

Let us consider now the classical equation of motion for a viscous fluid: 

dtipmUi) + dk^ik = 0, (6) 

where pm is the mass density and 11^^ the momentum flux density tensor [40] . 

Iljfc = p5ik + PmUiUk - p{dkUi + diUk), (7) 

where p is the shear viscosity. Taking the mean of this equation, observing that = in homogeneous 
and isotropic turbulence, we obtain 

2 

li^r, = IIj^j^ = = Pcff =V + gPtUrb, P = ^{t)p + K{t) , (8) 

where A(t) is the variable cosmological constant, a variable parameter of equation of state, and Poff 
is the effective pressure which takes into account that the thermodynamical pressure is augmented by a 
term (2/3)pturb associated with the turbulent energy density, 

Pturb = -pmUiU,; = -pmU^, (9) 

meaning that pturb designates a mean quantity. 

3 FRW equations using a fraction / of the energy as turbulent 
energy 

Let us now consider this fluid in a cosmological setting. As usual in turbulence theory we may 
start by decomposing the fluid velocity Ui into a mean component Ui and a fluctuating component u[, 
Ui = Ui + u[. On the other hand, in a comoving reference frame Ui = 0, so that we can simply replace u[ 
with Ui. Requiring only the turbulence part of the cosmic fluid, far from the non- viscous (non-turbulence) 
part, one can from now adopt the notation for which the non- viscous terms are subscripted by zero. Thus 
Eq. ([SI is rewritten as 

2 

Poff =P0 + gPtmb- (10) 

The turbulence energy in the comoving frame can be view as a nonrelativistic quantity, and the total 
energy density can be written as 

2 

P = PO + -^Pturb ■ (11) 
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Let us look for the cosmological equations when a constant fraction / of the energy exists in form of 
turbulence energy, and assume that the energy at initial time tin can be written as 

piUn) = PO (1 + /) . (12) 

From this, we see that / can be interpreted as the ratio between the turbulence energy and the total 
energy at Un, 

/ = (13) 
Po 

Since the rest mass is completely incorporated in the total energy, the condition / << 1 is required |24] . 
Therefore, for t > tin it is necessary pturb decaying with time as pturb according to the time 

evolution of the ordinary turbulence |24j . Thus, for t > tin the energy density can be written in the 
following form 

pit)=poit)[l + fp^it)] , p,{t)=(^^y . (14) 

The choice of pi (t) seems quick suitable since for t — tin , the relation (|12p is recovered. 
Let us now write the equation of state for the cosmic fluid as 

Po{t)^ujo{t)po{t)+Ao{t) . (15) 

In the same way as (jl2p . one can write the effective pressure as 

=w(Op(0+A(t), w(t) -cjoW[l + MW], A(t)=Ao(0[l + /Ai(t)] . (16) 

Through the condition / << 1, one can rewrite the effective pressure as 

Peffit) = po{t) + {woWpoW M) + Pi{t)] + Ao(t)Ai(i)}/ . (17) 

Analogously, the scale factor and the Hubble parameter can be expanded as 

ait) = ao(i)[l + /aiW] , (18) 
Hit) = Hoit)[l + fHiit)] . (19) 

The correction terms ujiit), Ai(t), ai(t) and -ffi(t) are all for zero order in /. Using (flQl) and deriving 
([TSl) one time with respect to t, one gets 

aiit) ^ Hoit)Hiit) . (20) 

From the first equation of Friedmann, = p, where we set k = SvrG = 1, G being the Newtonian 
gravitational constant, one gets 

Hiit) = ^piit) . (21) 



Hence, the Hubble parameter reads 

H{t) = Ho{t) 



6/5' 



(22) 



Now, we also have to determine HqII), and in order to do this, we need to use the equation of continuity 
related to the non-viscous components as 



po + 3i7o {pa +po) = 



(23) 



Let us now consider for simplicity the case where both ujq and Ao{t) are linear time dependent, i.e. 
u!o{t) = at + (3 and Ao(t) ~ jt + X. Therefore, the equation becomes 



whose general solution reads 
Ho{t) 



AHoiU 



B 



gC(t+Z5)^erf(t) - eC(*-+^)'erf(t„0 



where the input constant are defined as 



4 7 aV3 



B 



1 /27rV3 



Therefore, the Hubble parameter (j22p becomes 



B 



3C(t+i?)2erf(t) - eC(*-+^)'erf(t„) 



A 



A 



aC 
1 + 



1 . ft 



(24) 



(25) 



^+2nt 



6/5' 



(26) 



and the corresponding energy density reads 

3A2i?2(i»n) 



Pit) 



|B[eC(*+^)'erf(i) - eC(t.^+^)^erf(t„o] + a} 



1 + 1 



t 



6/5' 



(27) 



Observe that in general, for small tin/t (or t >> tin), the influence from the turbulence fades away for 
both the Hubble parameter and the energy density. 

Let us now analyse some special cases where cosmological feature may be obtained. 
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3.1 Analysing cuq as linear function of time and Aq a constant 



In this case, we have uJo{t) — at + /3 and Aq — A, meaning that 7 = 0. Then, by solving ([M)) . the 
Hubble parameter and the energy density are found as 



Hit) = 
Pit) - 
Z{t) = 



VXHoit^n)Z{t) 



^9aH'^it,n)it - t„) + XZ^{t,n) 
i\H^{t,^)Z\t) 



%aHl{t,n)it~t,^) + \Z^{t,n) 



t 

6/5 



6/5' 



(28) 
(29) 



Z_ 



1 /A 
3 



(ai + /3 + 1)2 



L J ^ {at + 13 + 1)'' 



where Z^ = CiJi, + C2Y^ is the general solution of Bessel's differential equation, Ci and C2 being arbitrary 
constants. One can observe that for t = — (/3 + l)/a both the Hubble parameter and the energy density 
diverge simultaneously, corresponding to the singularity of type-I (the Big Rip) . 
• Linear uJo{t) and vanishing cosmological constant 

In this situation the cosmological constant is null, Aq — 0, whilst uJo{t) — at + (3. Thus, by solving 
one gets the Hubble parameter and the energy density 



H{t) = 
Pit) = 3 



Ho{tin) 



Hoit^n) it - t™) [1(1 +f3) + ^it + t„) 
H^itin) 



1 



1 + 



{Ho{t,n) it - t„,) [|(l+/3)+ 3^(t + i„,)] +l| 



1 + / 



t 



6/5' 



6/5' 



(30) 
(31) 



These results generalize the ones found in 24 , which are recovered for a — 0. Like in vanishing turbulence 
case, finite time future singularity occurs at ts such that 



ai^ + 2{l + I3)ts-2{1 + t,n - ■„ + 



3Hoitin) 







For a 7^ 0, one gets 



± 



1 + /3 



Observe that, for a = 0, the singularity time readi 



tin 



3aHo{tt 



(32) 



(33) 



(34) 



3\1 + P\Hoit,n) ■ 

This later is exactly the result found for the singularity time in j24j . where they used 7 at the place of 
1 + (3 and (constant) at the place of f3. Consequently all the results found in that paper should also 
be found here. Since the goal of this paper is to search for more general cases, let us reconsider the case 



^This value has to be directly used in equation II32I I 



7 



where ujQ{t) is hnear time dependent and the cosmological constant vanishes. One can now search for the 
correction term ai which is such that di = HqPi/2. One gets 



aiit) 



1 



'2t.„y/' I 6/5... p 6 11 



-/2^'W2Fi (^^,^,^,g,{t)^ I ,(35) 



where the functions fi{t), and gi{t) are defined by 

At 



h{t) = 



B - y/B^ - 4AC + 2 At 



f2{t) = 



At 



B + \/B^ - AAC + 2At 



B - - 4AC 
B - VB^ - AAC + 2 At 



, 92{t) = 



A^—^ i?=-(l + /3), 



B + - AAC 
B + /B2~"4]4C + 2At ' 

■ - ^"4. - ^ (1 + /3) , 



(36) 



whilst the scale factor ao(t) in the vanishing turbulence case reads 



ao{t) = 



: arctan 



/ B + 2At 



(37) 



yjAAC - S2 V\/4AC- 
Since we are leading with the case where the cosmological constant vanishes, the correction term Ai is 
also null, (Ai(t) = 0). We can now proceed to the calculation of the correction term 0Ji{t). To this end, 
let us make use of the second generalized equation of Friedmann 



iH\t)-2H{t)=p,jf 



(38) 



The left hand side of ([55]) can be calculated from p^ . and comparing the result with the right hand side 
of (fT7|) . one gets 

1 



wi(t) 



Ho{t) + -SH^it) pi{t) + Hoit)pi{t) } , 



wo(i)/Oo(0 

{16A - 15a)t2 + (B - 15 - 15/3)t + 6C (t 



6/5 



\U{at + ;3) V * y ' ^^^^ 

Therefore, with an inhomogeneous equation of state of the form (jl6p . with vanishing cosmological con- 
stant, the correction term cji(t) associated to wo(t) is given by p9p. Observe that at singularity time, 
ts > tin, this correction term becomes inefficient such that the turbulence effect in uj{t) falls down and 
this latter reduces to i^o(0- 



3.2 Analysing Aq as linear function of time and uq a constant 

In this case we consider a = such that ujq = /3 and Ao{t) ~ jt + X. Thereby, the general solution of 
((24)) becomes 

6/5" 



3A 



2-' \ 4^-3A 
3 



(40) 



34"^ - 24 + (T 4 37 
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Here, it is important to note that this solution makes sense only if /? ^ —1, otherwise, we fall in the case 
where the Hubble parameter is a polynomial function of the cosmic time. But we will see this later. So, 
in this case of constant wq E^nd linear Aq, the energy density reads 

6/5" 



1 + / 



4e-3A ; 



(41) 



Since the parameter /3 is non-null, one can easily observe that both energy density and Hubble parameter 
diverge only for complex value of ^, consequently, for complex value of the cosmic time. Therefore, one 
can conclude that this case is interesting since the big rip cannot appear for finite real time. Note also that 
since the wq is constant, its correction term wi vanishes. Thus, one has now to search for the correction 
term Ai(i) corresponding to Ao(t). As we have performed in the previous case, now looking for Ai(t), 
one gets 

Ai(i) = {c^o(i)Po(i)PiW + 6i?2(i)i7^(t) +ijo(i)i/^(t) +i/o(t)iy^(t)| , (42) 

which, in terms of ^, yields 

Note that the sub-case where the parameter luq vanishes whilst Ao(t) is still linear function of time can 
be easily found from this later by setting /3 = 0. Is important to recall that the results are physically 
the same, because this is just a consequence of rescaling the parameter of Hubble by a multiplicative 
constant (leaving /3 + 1 to 1) in the equation (p4|) . Other interesting case to be mentioned is that in 
which both parameters ujq and A are constants. In this case, one has a = ^ = 0, luq = (3 and Aq = A. 
Expressions of the Hubble parameter and the energy density may be obtained by making use of the above 
assumptions in piI)) - (HT]) . In this subsection, one can conclude that for a constant and linear time 
dependent cosmological constant, future singularity can never occur, at least for real cosmic time. 

3.3 Treating uq as an oscillating function, with a vanishing cosmological con- 
stant 

Here, instead of using uJo{t) as linear function of time, we assume that it is an oscillating function of 
time. Let us consider ujoit) = — 1 + cos(fei). By making use of this later and vanishing the cosmological 
constant, one gets from the equation 
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and the Hubble parameter within turbulence reads 

bHo{tin) 



H{t) = 



h + iloHoiUn) COS [h{t + Un)/2] sin [h{t - U^)/2] 



where li is an integration constant, and the corresponding energy density is 



Pit) = 



6/5' 



1 + / 



6/5' 



(45) 



(46) 



|6 + iloHoiUn) cos [bit + Un)/2] sin [b{t - t„0/2] y 

It is also obvious here that for t >> tin, the turbulence effect fades away both for the Hubble parameter 
and the energy density. In order to analyse very clearly the occurrence of finite time singularity, let 
us reconsider the expression (j44p . Note that if |^i| < Iq, the denominator can be zero, which implies a 
cosmological singularity. However, if |^i| > Iq, singularity cannot occur. It appears that both BiG Rip 
and Little Rip may appear in dark energy within turbulence consideration. This kind of results has been 
found recently by Brevik et al |43) but not within oscillating inhomogeneous EoS. 

The correction term uji{t) corresponding to ujQ(t), can be found in the same way as the previous cases 

by 

2 



wi(t) 



ujoit)poit) 
1 

1 — Iq cos{bt) 



3H^it)Hiit) + Ho{t)Hi{t) + Ho{t)Hi{t) - pi{t) 



2 

5t 



■ cos{bti, 



6/5 



(47) 



which fades away for t » ti. 



4 Turbulence as additive dark energy and stability analysis 

In this section, instead of using the perturbative approach, we assume the total dark energy density 
Ptd as the sum of the usual dark energy pd and the turbulence energy pmrb 

Ptd ^ Pd + Pturb , (48) 

and the first equation of Friedmann is written as 

= Pd+ Pturb + Prad , (49) 

where prad characterises the energy density of the radiation. Treating the turbulence, the dark energy 
and the radiation as three component fluid of the universe, the corresponding continuity equations read 

Pd + 5H{pd+Pd) = Qi , (50) 

Prad + SH {prad + Prad) = Q2 , (51) 
Pturb + 3i/ {pturb + Pturb) = Qs , (52) 
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such that Qi ^- Q2 + Qz — 0. We also define dimensionless density parameters as 

Pd Prad Pturh ( 



The continuity equations in dimensionless variables reduce to 

= ix[x{l+LOd)+y{l+UJrad)+z{l+UJturb)]-?'X{l+UJd) + ^ , 
^ = 3y[x{l+LUd)+y{l+0Jrad)+z{l+UJturb)]-iy{l+^rad) + ^ , (54) 
^ = 3z[x{l+0Jd) +yi^+i^rad) + Z {I + UJturb)] -'iz{l+ UJturb) + , 

where iV = Ina is the so-called e-folding parameter. In general the interaction terms Qi, i = 1, 2, 3, are 
functions of the Hubble parameter energy density pj. One obtains the critical points by equating the 
equations in (|54p to zero. Next, we will perturb the equations up to first order around the critical points 
and analyse their stability. Only the stable critical point will be taken into account, i.e, the critical points 
for which all the eigenvalues of Jacobian matrix are negative. These points are interpreted as attractor 
solutions of the dynamical system. 

In this paper we will be restricted to pure physical situation. Because of the turbulence, the kinetic 
energy changes into heat and the heat then becomes radiation. Therefore, the interaction essential occurs 
between the turbulence and the radiation. Consequently, we may set the interaction term Qi to zero. As 
well known, we can make use of the value of EoS of radiation {ujrad ~ 1/3). Since the total content of the 
universe has to conserve, one gets Q2 — —Q3 = Q- In such a situation, since dark energy is assumed to 
not interact with the radiation, or the turbulence, one can consider its related parameter of equation of 
motion uJd as constant and more precisely about —1. Therefore, the dynamics of the system essentially 
depends on the turbulence and the radiation. Hence, the system ([54]) reduces to 



dN ^ 3H^ ' 



(55) 



dN — ""^^ (1 + ^turb) - 3773- 

We assume the interaction Q between the turbulence and the radiation as a quantity proportional to 
PradPturb/ H , i.c, Q — qpradPturb / , whcrc g is a real parameter. Thus, the system ([55]) becomes 

^ = -^V + <iy^, (56) 

jjf = -3z (1 + i^turb) - qyz ■ 

As inhomogeneous equation of state, we assume ujturb = kiH^^ pturb/i + ^2, i-e, ujturb ~ kiz + ^2, where 
ki and ^2 are two negative real constants. This choice of the sign of ki and ^2 is purely cosmological 
since as the turbulence contributes as a type of dark fluid, its pressure has to be negative in order 
to insure the acceleration of the universe and a necessary condition for this is assuming ki and ^2 as 
negative constants. The system presents two realistic critical points: a){yc, Zc) = [0, —(1 + fc2)/fci] and b) 
[— 3(q + qk2 + 4fci)/q^, 4/g] . The first critical point correspond to the situation where the kinetic energy 
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of the turbulence is still conserved, i.e, there is any change of the kinetic energy into heat. However, the 
second critical point is related to the situation where there is a real dynamics of the system. In this case, 
part of kinetic energy is lost, yielding radiation. We will focus our attention to this second critical point 
and analyse the stability of the system about it. To this end, we will consider small perturbations around 
this critical point and write any point (y, z) of the system as {y + Sy, z + Sz). Thereby, the dynamical 
system becomes 



By denoting M the matrix associated to the above system, one can find the eigenvalues by solving the 
equation det [Ai — Xl\ =0. The solutions of this equations read 



It can easily been seen that the sign of the eigenvalues strongly depend on that of the interaction parameter 
q. One can observe two general cases: g < and q > 0. For q < and k2 < —1 — 4fci/g, one gets 
—6ki/q < and both Ai and A2 are negative. In this case, the critical point is an attractor and the 
system may be stable. Still for (7 < 0, if —1 — 4A:i/q < fc2 < 0, Ai > and A2 < and the critical point is 
called saddle point and the system cannot be stable. For q > 0, —Qki/q > 0, always one gets A1A2 < 0, 
or Ai > and A2 > 0. In such situations, the system cannot be stable. 

We clearly see that the stability of the system is realized only when ki/q > and ^2 < — 1 — 4fci/g. By 
the use of these conditions, we from (|57p that as the universe expands, i.e, for increasing N (increasing 
scale factor), the perturbation of the turbulence and radiation energy densities decreases. 

5 Conclusion 

In this paper, we investigate the effects of turbulence in dark fluid universe at late accelerated ex- 
panding universe. We focussed our attention to the model of linear inhomogeneous EoS, with linear and 
oscillating dependence on time. We included turbulence effect by the addition of a constant fraction to 
the laminar ordinary energy density. We then performed a perturbative analyse checking the possible 
appearance of future singularity. The obtained results generalize that of Brevik et al [24. Moreover, for 
each considered expression for the parameters and A(t), we determined their respective correction 
terms. 

In the second step, instead of considering the perturbative approach, we considered the dark energy 
contribution as a sum of a laminar part of dark energy and the turbulence part. Then, we study the 
stability of the system, considering that the dynamics of the system essentially depends on the turbulence 




(57) 




(58) 
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and radiation energy densities. The results shows that the stabihty of the system is obtained when the 
turbulence energy density decreases as the universe expands, while radiation energy density of radiation 
grows. Therefore, when turbulence contribution is taking into account in the total dark fluid, according 
to the statement that kinetic energy changes into heat, the stability is always realized. Hence, precisely 
in phantom region, due to this probable stability, one may conclude that turbulence may avoid future 
finite time singularities, and this when all the input parameters g, ki and /c2 are negative and satisfy the 
relations fci/g > and fc2 < — 1 — 4fci/g. 

Acknovifledgement: M. J. S. Houndjo thanks CNPq-FAPES for financial support. M. E. Rodrigues 
thanks UFES for the hospitality during the elaboration of this work. C. Ai'namon would like to thank 
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